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SIDESLIP OF WING-BODY COMBINATIONS 


By Paul E. Rubbert 


THE BOEING COMPANY 
Seattle, Washington 


1.0 SUMMARY 


A theory is developed for predicting the aerodynamic properties of an airplane in side- 
slip. It is developed around a model problem consisting of a wing-body combination. 

The basic approach used is perturbation theory. A solution is assumed in the form of an 
asymptotic expansion in powers of the small parameters governing angle of attack, angle of 
sideslip, wing camber, wing thickness, body radius variation, and body camber. Second-order 
terms governing the interaction between sideslip and the other parameters are retained. It is 
found that when wing dihedral is large, the dominant sideslip effects occur in the first-order 
terms. When dihedral is small, however, the dominant sideslip effects occur at second order 
and depend on products of sideslip and angle of attack, sideslip and wing camber, etc. 

Particular integrals for the inhomogeneous second-order problems are given which auto- 
matically satisfy the wake boundary condition. Hence, the remaining homogeneous part can 
be formulated in terms of distributions of sources and elementary horseshoe vortices on the 
boundary surfaces in a manner entirely analogous to that commonly done for the first-order 
problems. 

The present formulation is shown to produce results which reduce to known solutions 
for the infinite yawed wing and for a general planar wing in sideslip. 
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2.0 INTRODUCTION 


The development of a small-disturbance theory for the aerodynamics of an airplane in 
sideslip entails a description of the airplane geometric features and the inclination of the on- 
coming stream in terms of small parameters. The parametric description must be such that 
the flow field reduces to a uniform stream in the limit as all parameters vanish. 

For planar wings the problem can be formulated using coordinate systems aligned either 
with the freestream or along the plane of symmetry of the wing as shown in figure 1. 


£ 



(a) Skewed Planform 




X 

(b> Skewed Freestream 


FIGURE 1. -ALTERNATE FORMULATONS 


In the former case the planform is skewed with respect to the coordinate system by the 
sideslip angle & , which remains constant in the limit as the angle ol attack and wing thickness 
parameters approach zero. Thus, /£ does not appear directly in the. problem as a small param- 
eter and is not subject to a limit process. However, in the case of the skewed stream formula- 
tion,^ appears directly as a small parameter for which a limit processs may be defined. 

Lowest order solutions from the skewed planform formulation display the dominant 
effects of sideslip. However, since sideslip does not appear explicitly as a separate small 
parameter, the rates of change of flow properties with angle of sideslip, and lateral stability 
derivatives, must be computed by finite difference from solutions obtained at different angles 
of skew. 

The lowest order solutions in the skewed freestream formulation do not produce the 
dominant sideslip effects. The first order in^ solution for a flat wing with no dihedral pro- 
duces nothing, since the sideslip component of the freestream, which is directed along the y 
axis, remains undisturbed by the presence of the planar wing. The dominant sideslip effects 
occur in second-order terms involving products of sideslip and angle of attack, sideslip and 
thickness, etc. These terms contain the sideslip parameter explicitly, and provide directly the 
rates of change of flow properties with sideslip and the interaction between sideslip and angle 
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of attack, etc. Thus, the results produced by this formulation appear in a form that is more 
directly applicable for stability and control work and other aerodynamic applications. 

For more complex configurations containing a fuselage, nacelles, dihedral, etc., a formu- 
lation using the skewed planform approach becomes very unwieldly, because of the difficul- 
ties involved in defining parametrically the geometrical shape such that the flow reduces to a 
uniform stream in the limit as the small parameters vanish. The skewed stream approach, on 
the other hand, can use the same parametric geometry definition normally used for symmet- 
ric conditions. Hence, for the present work a decision was made to use the latter formulation, 
even though it entails a higher order development. 
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3.0 SYMBOLS 


AqO 

aR(x) 

B 

s 

C Pa 

Cp a* 

C P C 

C Pc^ 

C P* 


V A/S 

% 

C P^ 

c p e/3 

c p^ 
C P 1-6 

cG(x) 

K s 

K V 

JL 

M 

m 


freestream sonic speed 

body radius variation from the mean body surface (eq. 1 ) 

\! | 1 - M 2 | 

complete pressure coefficient (eq. 1 0) 

pressure coefficient due to body radius variation (eq. 28) 

pressure coefficient due to body radius variation and sideslip interaction (eq. 28) 
pressure coefficient due to body camber (eq. 28) 

pressure coefficient due to body camber and sideslip interation (eq. 28) 
pressure coefficient due to angle of attack (eq. 28) 

pressure coefficient due to angle of attack and sideslip interaction (eq. 28) 
pressure coefficient due to sideslip (eq. 28) 
pressure coefficient due to wing camber (eq. 28) 

pressure coefficient due to wing camber and sideslip interation (eq. 28) 
pressure coefficient due to wing thickness (eq. 28) 

pressure coefficient due to wing tliickness and sideslip interaction (eq. 28) 

body camber distribution (eq. 2) 

potential induced by a point source (eq. 36) 

potential induced by an elementary horseshoe vortex (eq. 34) 

lateral coordinate along a surface 

freestream Mach number 

source density 
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n 

rT 

P 

q 

r 

r 

S 

s 

u, v, w 

Voo 


x, y, z 


OC 

t 

AQ n 
A q t 


£ 

0H(x,s) 


Jl 


/> 


coordinate normal to a surface 
unit vector normal to a surface 
radius of mean body surface (eq. 1) 
local flow speed 
body radius 

radius vector from a surface point to a field point 

denotes a surface 

spanwise coordinate 

perturbation velocity components 

freestream speed 

Cartesian coordinates (fig. 2) 

angle of attack 

angle of sideslip 

vorticity density (also the ratio of specific heats) 

discontinuity in the normal velocity component across a surface (fig. 1 1 ) 

discontinuity in the tangential velocity component in the lateral direction 
across a surface (fig. 1 1 ) 

lateral width of a surface strip (fig. 1 2) 

wing camber form (eq. 3) 

sweepback angle 

see figure 6 

local coordinates (figs. 12, 13, 14) 
density 
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r F(x,s) 


wing thickness form (eq. 3) 
complete velocity potential 
<f p denotes a particular solution 

y,s portion of a perturbation potential contributed by sources 

y,v portion of a perturbation potential contributed by vortices 

f wing dihedral angle (fig. 4) 

gradient operator 

[ ] denotes discontinuity of a function across a surface 

Perturbation velocity potentials (eq. 5) 


% 

% 

% 

4 

% 

% 

% 

%o 


angle-of-attack effect, 
first order sideslip effect 

axial flow over uncambered body with a thin, uncambered wing. 

body camber effect. 

wing camber effect 

wing thickness effect. 

angle-of-attack and sideslip interaction. 

body thickness and sideslip interaction. 

body camber and sideslip interaction 

wing camber and sideslip interaction. 

wing thickness and sideslip interaction. 
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4.0 GEOMETRY DEFINITION 


The basic theory is developed around a model problem consisting of a wing with arbi- 
trary dihedral attached to an infinitely long fuselage. The results can be readily extended to 
include multiple lifting surfaces such as horizontal and vertical tails, as well as nacelles, struts, 
and similar appendages. The effects of body truncation in the form of a pointed nose and tail 
can be added in a manner similar to that described by Woodward (ref. 1), wherein the classical 
slender body solution representing the pointed nose and tail is patched into the results of the 
present analysis. 

Let all length dimensions be scaled by the average wing chord and assigned to be of order 
unity. The X axis of the coordinate system is aligned lengthwise with the configuration, as 
shown in figure 2. 



FIGURE 2. -COORD! NA TE SYSTEM 


The y axis points to the right of an observer seated in the body, and the 3 axis is directed 
upward. 

The freestream velocity vL is inclined at an angle ac to the X ? plane. The component 
of vL in the * y plane is inclined at the sideslip angle j6 to the X axis. When the compo- 
nent of ^ along the y axis is directed in the positive ^ direction, shall be defined as 
positive. 
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The body shall be of circular cross section, and may be cambered. The body radius may 
vary slightly along its length. The nomenclature chosen for the body is shown in figure 3. 

H 



The radius distribution of the body is denoted as 

r ~ P + a. x) (1) 


where 

P = constant =0(1) 

R (X) = O CO 

a. << I 

The body camber distribution is 

2. = C G Cx) ( 2 ) 


where 

G OO - o< \ ) 
c <.< i 

With this notation, the body approaches a cylinder in the limit as a., C - — O. 

The wing shall be positioned along a mean surface described by Z = Z W (S), 
(S ) } whose generators are parallel to the X axis, as shown in figure 4. 
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Rear View 

FIGURE 4 -WING DIHEDRAL SHAPE 

The wing dihedral distribution is defined by tf> (&) , where 5 is distance measured along 
the span. For the initial analysis, ip fsj is assumed to be 0(1). 

The wing surface shape is described parametrically as a function of the coordinates 
(Xj 5} . Decomposed into camber and thickness, the wing surface shape becomes 

y(Xj 3 ) » (s) + F (x, y - elUXj s) sin p 

(3) 

1(X3) = z w (5)t TF(k ; 5) Cos 9 + OHCXjS) cosy 


where 

RXj s) describes the thickness shape 

U ( x ? 3 ) describes the camber shape 

& are small parameters governing the magnitude of the thickness and camber 
distributions 

The upper and lower signs refer to the upper and lower wing surface, respectively. With these 
definitions, an airfoil section generated by a cutting plane normal to the mean wing surface 
assumes its usual definition of 


^2-0 ~ ~ 7* F ( X f S) + &H ( XjS) 

AW Foil 
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5.0 EXPANSION OF THE VELOCITY POTENTIAL 


An asymptotic expansion of the velocity potential containing all first-order terms plus 
the second-order interaction terms involving the sideslip angle & is assumed to be of the form 

i = + <* V"o +6 V, + * % + c % 

■he W * r ?s + *6 + * a y 7 + Ctf y s (5) 

+ y B + Tjg y to ) -t oca > stc -) 

The various perturbation potentials appearing in this expression can be interpreted in the 
following manner. 


% 

% 

% 

% 

% 

*4 


angle-of-attack effect, 
first order sideslip effect 

axial flow over uncambered body with a thin, uncambered wing. 

body camber effect. 

wing camber effect 

wing thickness effect. 

angle-of-attack and sideslip interaction. 




body thickness and sideslip interaction. 


% 


body camber and sideslip interaction 
wing camber and sideslip interaction, 
wing thickness and sideslip interaction. 


The objective in choosing this particular form is to obtain a sequence of simplified prob- 
lems governing the primary effects of camber, thickness, angle of attack, and sideslip. Second- 
order terms involving the sideslip angle /& are retained in recognition of the fact that the 
dominant sideslip effect may appear as either a first-order or a second-order effect, depending 
primarily on the magnitude of the wing dihedral. For the case of a wing with large dihedral, 
the potential % will be of order unity (its magnitude is governed primarily by the wing 
dihedral). In that case, the first-order term /$ tfj gives the dominant sideslip effect and all 
second-order terms could be neglected. For the case of a wing with no dihedral, however, 
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is zero and the dominant sideslip effect is given by the second-order terms. Thus, in the 
general case it is necessary to include at least some ot the second-order terms given in eq. (5) 
to ensure that the dominant sideslip effect is not overlooked. Arguments will be given later 
for eliminating some of these terms. 

Substitution of eq. (5) into the equations of motion and the boundary conditions will 
yield, after equating like powers of the small parameters, a series of equations and boundary 
conditions for the determination of each of the separate perturbation potentials. 
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6.0 BOUNDARY CONDITIONS 


The boundary conditions consist of the requirement that the flow be parallel to all solid 
surfaces, the Kutta condition for subsonic trailing edges, continuity of pressure across the 
wake, and the vanishing of disturbances at infinity. The requirement that DB/Dt = 0, where 
2 i) defines the position of the solid surfaces, leads to the appropriate expressions 
for the solid surface boundary condition. 


6.1 BOUNDARY CONDITIONS ON THE WING 

The equation of the wing surface, (3), together with eq. (5) for the velocity potential, 
are substituted into the expression DB/Dt = 0 to yield the result 

[-r-K + [»+*<&* + ^vix + + ?,oA 

4- [ •+• 5 1 n ^ — 'X (-§§- cos <p -F ^ris SLr i y) +■ e ^s"s cos ^ 

- n sLti )][/5+ oc V> oy 4/5 ^ 2 y + ‘ + f lo y ] ( 6 ) 

+ [±CM ^ - r Si n V 4 F Co -s <P + 6 ( fs~ 5Ln V 

4 M cos «p)] [of -Erf 4 /3 f /a 4 <2 +• • ‘ + T ^ tfo, ] = ° 

0*6 = 0 

where higher order terms not considered in the expansion (5) have been deleted. 

The various potentials on B - O must now be expressed in Taylor series about the 
mean wing surface. Since we are interested in retaining only those second-order terms 
involving products of , this operation will introduce extra terms only trom $ 
appearing in eq. (6). The proper expansion is 

±e» r f( KS) (7) 

+ o(~t * e 2 } T q) 

and 

(b = ±o) = ^,2 r?7= ?y, vr~ o; ^ T r< " Y/ s; 

4 o (r 2 , ©*, , 

where Yl is a normal coordinate measured upward from the mean wing surface. The ± sign 
denotes the upper or lower surface of the wing, respectively. 
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Inserting eq. (7) into (6) and equating orders of magnitude finally results in the boundary 
conditions to be applied on the mean wing surface for the various potentials. 


o(a): = 

377 

- Co 3 <j> 

o (<j): s = 

3>?? 

S/* ^ 

O (cl) » ^ ^ 

= o 

o (a) : #-&- 
' S> n 

» o 

o (e) : s ^ 
9 n 

o>W 

a x 

o(r),y% 

+ 5F_ 
S X 

ft ft 

1 

0 

- o 

o <•<*; : If 

= o 


- o 




+ 


3Ai 

3X 


+ Hsiy? tfj *-VyL - Vcos V % 

dn 9 71 


° ^ T/3) : ITrT “ "*■ * A - ^ c<?s v) + ** /2 ( F sin V) 


+ <f 


/x 


o> F 

d y 


■4- F si?? ^ 

3 77 


- r cos <*> 


g **/* ? 

s>* / 


where 


■3 71 


Cos 



St 7j 
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Note that the potentials appearing on the right-hand side of the and 

expressions may be discontinuous across the mean wing surface, indicating that a source-like 
term may be present in the solution. 


6.2 BOUNDARY CONDITIONS ON THE WING WAKE 

The condition to be satisfied is that the pressure be continuous across the wake surface. 
Introducing expansion (5) into the exact expression 

c p"£r* <? -*.')]$> -ij (9) 

where 

M = freestream Mach number 
V - freestream speed 

OO 

^ = local flow speed 

A ^ = sound speed in the undisturbed stream 
and expanding in powers of the small parameters yields the following expression for Cp 

Cp = -2 ^ V»„ t /S <f, x + * f zx + c V> 3 X * s * r tis* 

+ at/cr^ + r ajf + y, z + v % ■ v f, - m* ? 0K </>„} 

+ a 0 (<fi rx + if tj/ x- V f, • v V,* fix) 

+ C 4 ( t»x 1- * V Y, ■ V V 3 - f,x %X ) <l0) 

* + fig I- vy ■ fix'?**) 

+ ( % ox + Vs 3 + • 7 A - " 2 f, x 9, x )J 

The pressure coefficient on the actual wake surface (whose position is generally unknown) 
is next expressed in a Taylor-series expansion about the mean wake surface, a cylindrical sur- 
face lying directly behind the mean wing surface. Let the deviation of a point on the actual 
wake from the mean wake surface be denoted as Ay (* ) ^ A*Zr (k ) f as shown in figure 5. 
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The distances A y (x) ; A £ (x) are functions of the small parameters and may 
be expressed as 



+ ^ Z + 

CA^ 

+ + T A *5 

-h 

0 ( OL/S ) J 

e-t c.. 

(ID 

AZ (x) - oi A2 o +^AZ, 

-h CL A 2 Z 

+ CAZj 

4- -h 

+ 


c. 



With this notation, the values of the potentials or their derivatives on the actual wake surface, 
expressed as a series expansion about the mean wake surface, become of the form 


tL 


dc 

w a ice 


~ /mean 

Wake 


2><P / Au / 

£><y ' * A # / * 


A? -1-0 


mean 


fptCLTf 


Ayiij (AEfy 1 


( 12 ) 


Equations (1 1) and (12) substituted into (10) furnish the appropriate expression for the 
pressure coefficient on the wake. Imposing the requirement that the pressure coefficient 
must be continuous across the wake to all orders of magnitude then results in the following 
sequence of boundary conditions to be applied on the mean wake surface. (The notation [ 1 

is used to denote the discontinuity of the function enclosed in brackets across the surface.) 
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0 (<* ) •• [ f OK J - ° 


0(4): =0 

°<“> ■ [fix ] = O 

o(c) : J = o 

•’ [%x ] ‘ ° 

o (T) : [ ?sx J » o 

O (<X0) : [ ] = M * f /x J ] " [f/B ] ' ^ f ] 

A it° [V/Xff ]"*<?' [Vox* ]~ Ai °[fixe ] 

' * z < o. « 2 ] 

oU6), K ]= ^ 0„ •»** ]-fa,]- bn • rr.l (i3) 

A it/ [flX& ] - A tf [fgXi] 

0(c/SJ : [*> ax ] = M 2 f«, <*, X J -[f 3 s]-[ v f - 7 f 3 ] 

A V/[f,xy] ■ AZ -l[flXiJ 

°(o& ■ [fix 1 = m2 [r<* f*y] - [«*»] - f v f ■ 17 M 

- - a, K» ] - Ai ' [r„ E ] 

o( Tx3) :[ ¥> o J = M ? [V, X <fi u ] - [<e s j - [ V f, ■ V ? 5 J 

- A j/' [fsxy J - A *'[fsxe ] 
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At first glance it appears that it may be necessary to compute the wake position from 
the first-order solution in order to evaluate , and which 

appear in the second-order problems. However, it can be shown that the first-order X -velocity 
components and their derivatives are continuous across the 

mean wake surface, and hence the terms multiplying 4 , dy, , Ql a ; and 

are all zero. The second-order wake boundary conditions are, in fact, independent of the 
wake deformation. 


6.3 BOUNDARY CONDITIONS ON THE BODY 
The equation of the body surface is 

3 (X; y, Z ; = 0= Z - CG (x) - y 2 (14) 

where 

rex) = P + <2 fx) 


The requirement that DB/Dt = 0, necessary to produce the solid surface boundary con- 
dition, leads with the use of eq. (5) to the expression 

- [ c 7T + ^a 7t][' + * ‘ ' ' + r> " 5 " J 

-+ [/&+ * <P 0 U +0 Vtg * • ' ' + ^ fto# J 

(15) 

■h [a + dL +£ +■ * • • + 't'fi Koz] = G 

o* y , z) = o 

, as sketched in figure 6. 

The next step is to express the velocity components on 6 (x t y , z ) = O in terms of 
their values on the mean body surface r* P by means of a Taylor-series expansion about 
r * P . Figure 7 shows the spacial relationship between a point A on the actual fuselage 


where _ f y 

A Vr*- 
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FIGURE 7.— EXPANSION OF BODY BOUNDARY CONDITIONS 

surface S ( x, ^ , Z )- O and the corresponding point A’ on the mean fuselage surface ot 
radius P about the X axis. The velocity potential at the two points is related by the 
expression 

<P (A) - . [ C G (x) + 4. *( X) Sir?^ 1 

+ • a. R(x) cos/u + o { c 2 , a * (i6) 

s f(A') + *1"*? • a. R(k) +. CG Cx), 

9 ?t ~ 2 

where Ti is a coordinate normal to the mean body surface and positive when directed outward. 
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Each term in eq. (15) is expanded in this manner. Finally, terms of equal order of mag- 
nitude are equated to produce the body boundary conditions for the various potentials. For 
convenience, certain higher order terms can be retained in the first-order boundary conditions 
to facilitate the numerical application of the Karman-Moore method (ref. 2). In the following 
equation, which follows from eqs. (1 5) and (16), these higher order terms, which in principle 
are not needed but in practice may be convenient, are underlined. 


0(pC ) : 


a ¥><, 

g n 


JR , 

Stn jLi. + d x 


°(/3) : 


s> _ _ 

3 V 




0(4) : 

7 Z 


o(c) 


a ¥> a 
a n 


dR 
d X 

o/ G 
dX 


dR 

d X 


! 2K 


sin ju 


4 t* 
d X 3* 


0(0) : g V* 
a * 


d x '4-x 


°( T ) •• - a gg- ts x (l7) 

^ (cCfi) : « o 

o> 77 

•• - *,* (fir -*<*'> t* cos * + **/<) 

o(cjs) r |fi - f, x jf- -&<*> 55 t>/y Cos ^ +K* at */‘J 

0(94) : = o 

3 77 


°( r* ) ■■ 

d 7? 


= O 


where — ^ 


COS/4- -h SiTfyO. 


on Vy 2 y- 2* — P =• o 
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7.0 FLOW EQUATION 


The equations of motion and continuity for an irrotational flow can be written as: 


V ({ i 2 ) + ^ VP''° 

V * (f V O 


( 18 ) 


The assumption of irrotationality is valid to the order of approximation sought. Eliminating 
the pressure from these equations gives 

2A 2 vV = (V^ Z )‘ (V?) (19) 

where 

i (S-!)( Vl - % z ) 

Substituting the expansion (5) into eq. (19) and equating orders of magnitude produces 
the governing equations in a straightforward manner. The results are given in the following 
section, which summarizes the equations and boundary conditions for the various potentials. 
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8.0 SUMMARY OF THE COMPLETE FIRST- AND SECOND-ORDER PROBLEMS 


The following separate problems are obtained by grouping together the results given by 
eqs. (8), (13), and (17) and the results of the previous section. The boundary conditions 
listed are to be applied on the mean wing surface, the mean wake surface, and the surface of 
the mean cylindrical body. The normal coordinate, n , is positive when directed upward from 
the mean wing surface or outward from the mean body surface. 

O(oc) ■ 


Differential equation: 0- M 7 ) ft 


+ U> 

o xx ' a 




fc 




Boundary 

conditions 


/Wing: 


'wake: 


\body: 


J5J1 

3 7 ? 


Cos <p - <fou V m “ Co5 * 


K] - ° 

^ " fou Cos M- 3*- ~ ^ d~x 




0(4): (>-mVv ixx + Vvy* 

wing: = Sen <P 

drj 

wake: \? fK ] = ^ 

body: =- c ° 5 ^ 


0 ( 0 .) : 


O(c) : 


(/-M 2 ) ft 
3 ft 

Wlne: T7 

XX + 

= o 


wake: 

] 


hody: ^ 

_ d£ 

Jx 

' ^ ^ " ft 

dx ’zx 

0- 

^3 XX 


s Pm 

wing: - J 

a 7? 

- o 


wake: 

] = ° 


bedy 

o/ <£r 

dx 

dR /j0 
V- «• J X 


( 20 ) 


21 



(Pi's) .• 

ra- 

^■^^4-xx ^ $ 4 . yy +~ 

(p = o 

r ^2Z 


wing: 

o >¥>4 _ S>// 

5> 7? 5> X 



wake: 

0 

It 

>< 

51. 



body: 

3 ^ ^ a Zt* f 

dr? c/ X 


<P^r; 
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Note: The underlined terms are of higher order and can be deleted if desired. 
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9.0 SIMPLIFICATIONS 


The complete expressions for the second-order terms summarized in the previous section 
can be simplified considerably without compromising the basic objective of determining the 
dominant sideslip term or terms. It was previously pointed out that different terms may domi- 
nate, depending on the particular aircraft geometry. Thus, for example, if the wing dihedral is 
large, the term will give the major sideslip effect, and it is unnecessary to take into con- 

sideration the higher order terms. On the other hand, if the aircraft geometry is such that the 
size of the term /S </>, is small (for example, /$ /, is zero for a planar isolated wing), then it 
is necessary to include the second-order terms. The clue that enables one to effect a simplifi- 
cation is to note that, for the example cited above wherein /3 <P, was zero, all of the terms 
involving vanish in the second-order problems, which simplifies them considerably. 

Let us first examine in detail the o(fi) and o(oifi) problems, which are reproduced 
below for convenience. 
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The magnitude of the first-order term ^ is governed by the boundary conditions of 
the o (g) problem. The wing contributes an effect proportional to , and the body an 
effect proportional to p z , the square of the mean body radius. Hence, the magnitude of 
the & •P, term may be denoted as o (fi ,/Sf >7 \ If and p 1 - are small (ot the same 
order as oc , say), then the second-order terms are comparable in magnitude to this one. 
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Turning now to the Ofa/S ) problem, we split the particular solution of the inhomo- 
geneous differential equation into two terms, and , satisfying the equations 

O-rtV f P , xx +<fp,„ + f P , a - fexy + 


O-m - 2M * w (? v»„ • vt* * m 2 )j x 

where 


^ particular " P t 

solution 



(24) 


An estimate of the magnitudes of Pp, and Vpz is obtained by examining the size 
of the inhomogeneous terms. The flow due to angle of attack, P 0 , will be 0( 1 ). The term 
P, is p 2 ) . Hence the magnitudes of the first-order term and the second-order 

particular solutions must be as follows: 

tv, = o(/g V) pp*) 

°<fi <fr, ° (<*/*) ( 25) 

<*/S = 0 (*4S 9>, 

The simplification that can be achieved is now evident, for note that the last term, 

Ci jS P p , is always an 0(oL) smaller than Pi , regardless of whether ^ and P are 

large or s’mall. Furthermore, when V and P z are small enough to be of the same magnitude 
as oi, , which is the condition for which aC/S Pp, must be retained, then the term *0 PP’z 
may be neglected in comparison with oifipp^ . Hence <*fl V *p z is always an oO*) 
smaller than the largest term, whether it be jS P, or <*0 P p, , and may always be neglected. 
One could also neglect that portion of Pp, arising from the term Pi on the right-hand 
side of eq. (23), since it too must always be an o(dL) smaller than the largest remaining term. 
However, it will be retained in the later analysis, since it lends a degree of symmetry to the 
relationship between oi and /S , 

Turning now to the boundary conditions and applying the same reasoning as above, it 
becomes apparent that the product terms involving Pt or its derivatives that appear in the 
boundary conditions for P* may also be neglected, for they too will contribute terms which 
are always an ofa) smaller than the largest remaining term. Hence the simplified problem for 
the potential governing the interaction between angle of attack and sidelip reduces to 
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wing: ^ - o 

wake: [£< J = " fag] " [«.] (26) 

body: | % -0 

where the underlined terms could be neglected if desired. 


The remaining second-order problems reduce in a similar manner, For example, consider 
the O(Gfi) problem involving , which is reproduced below for convenience. 
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Following the same argument as betore, it is apparent that the contribution to the par 
ticular solution of &A arising from the product terms involving on the right-hand side 
of the differential equation are always an o(e) smaller than the term fi <P ( , no matter whether 
and p*are large or small. When Vand P 2 are small enough to be of the same size as©, which 
is the condition for which the second-order term must be retained, then the product 

terms are an o(y , p* ) - a(e ) smaller than the remaining inhomogeneity involving 

c/> 4 * . Hence in all cases the product terms are an o(o) smaller than the largest remaining 

term and may be neglected. 


The same reasoning applies to all terms involving </>/ which appear in the boundary con- 
ditions. When ¥> and p 2 - are large, the entire second-order term is negligible, and when V 
and p 2 are the same size as & , then the terms in the boundary conditions involving *P, 
and its derivatives are an o(s) smaller than the remaining terms. It follows that all terms in 
the boundary conditions for which involve V* / or its derivatives can always be neglected. 
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10.0 SUMMARY OF THE REDUCED FIRST- AND SECOND-ORDER PROBLEMS 

In summary, the problems for the various potentials reduce to those listed below. The 
boundary conditions listed are to be applied on the mean wing surface, the mean wake sur- 
face, and the mean cylindrical body surface of radius P . The normal coordinate, n , is posi- 
tive when directed upward from the mean wing surface or outward from the mean body 
surface. The angle /4 is as shown in figure 8. The reduced differential equations and boundary 
conditions arise from application of the reductions derived in section 9.0 to the set of eq. ( 20 ). 
The respective expressions for the pressure coefficients come from eq. (10). 
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FIGURE 8.-CROSS SECTION COORD/NA TES 
Summary of Reduced Problems 
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The solution for the flow about a wing-body configuration in sideslip is given by the 
sum of the solutions yielding the individual perturbation potentials (c.f. eq. (5)) . 

It is interesting to attempt a physical interpretation of the second-order problems. Con- 
sider first the inhomogeneous terms appearing in the differential equations. They apparently 
serve to rotate the compressibility axis (the X -axis) through the angle fi to correct for the 
fact that in the first-order formulation the X -axis was not aligned with the freestream. This 
is apparent from the particular solution which for the o(ot0) problem may be written down 
by inspection as 


^ particular = * 0* " * fey < 29 > 

with similar solutions applicable for the other problems. The terms in eq. (29) are just the 
first-order terms in the Taylor-series expansion of required to construct the analytic 
continuation of the potential at any given point to a different point obtained by a 
rotation through an angle 0 . 

The wake boundary conditions of the second-order problems require the addition of a 
spanwise component of wake vorticity equal in strength to multiplied by the axial vorticity 
component generated by the first-order solutions. The result is to incline the vorticity vector 
in the wake to an angle /S with the X -axis, which is physically correct. 

The boundary conditions on the wing for the o^>nd *(t0) problems account for 
the effect of spanwise camber and thickness variations encountered by the sideslip component 
of the freestream. For many configurations the spanwise camber and thickness slopes are small 
enough to be neglected. 
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11.0 PARTICULAR SOLUTIONS OF THE REDUCED SECOND-ORDER PROBLEMS 


The distinguishing feature of the second-order problems that renders them more diffi- 
cult than the first-order problems is the presence of inhomogeneous terms in the differential 
equation. The first step in solving these problems is to establish the particular solution of the 
inhomogeneous equation. 

Since the inhomogeneous equations for the perturbation potentials % *f lo are all 
of the same basic form, one can search for a generalized form of the particular integral to be 
applied to all five cases. Let us describe the general problem in the following manner: 

(l-MV + <f>„ = 2M 2 (30) 

where 

ip = any of the second-order perturbation potentials 

cp ( °) = the corresponding potential appearing in the inhomogeneity 

All of the inhomogeneous equations appearing in (28) are of ttiis general form. The general 
wake boundary condition is written in the form 

[fx l-'-LfP’} on the mean <30 

wake surface 

One form of a particular solution of eq. (30) can be found by inspection to be 

= y - y (32) 

However, this is not the most convenient form, for it introduces additional singularities into 
the surface integrals appearing in ip^ ^ (which will be expressed in terms of source and vortex 
distributions on the boundary surfaces). A more suitable form for the particular solution is 
available, and will be described in section 1 1 .2. 

11.1 DISCUSSION OF HOMOGENEOUS SOLUTIONS 

As an aid in interpreting the behavior of the particular solution to be described in 
section 1 1 .2, it is convenient first to review the character and behavior of the homogeneous 
solutions of eq. (30), which are identical in form to the solutions of the first-order problems 
and whose behavior is well known. Homogeneous solutions will be expressed in the form of 
source and vortex singularity distributions on the mean boundary surfaces. Let us introduce a 
two-dimensional coordinate system X t £ on a mean boundary surface, where X ( is the 
axial coordinate and j£ the lateral coordinate, as shown in figure 9. 
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Mean boundary $urface,S 


7 : 



A third coordinate, n , is defined normal to the surface, having positive orientation when 
directed along (Xj x]T) • With this notation, the homogeneous solutions appear in the form 
of surface integrals, as given below. 


Vortex Distribution 

V (x> y, ?■) - *Yx,, Q.) k V (X J y, Zj £ )c/s ( 33 > 

where 

X , y- , Z = coordinates of a field point 

v 

( Xf'j.'i.) = portion of the perturbation potential contributed by vortices 

IT (y(, = vorticity density, on the surface, of strength established by the boundary 

conditions 

ds - d Ai 


The function X,, £) denotes the potential of an elementary horseshoe vortex 

located at ( *, y £) . It is defined as: 
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x-x,(s) 
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( 34 ) 
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(34 cont.) 
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h (rx-x,r5)) rz-2,rs)) 2 ) 

7(*-x,(s]) 2 - a 2 (T^TT^] 
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"r" = radius vector from (x, (s), (s), 2,(s)) to a field point (V, y, z) 

B =y /| l- M 1 1 

... , f 1 for U?-o| 

|o (lor co <L oj 


Source Distribution 

V^X,2,Z) = /^ £)K S (y, y, Z, x t> jl)Js (35) 

where 

</X, y, Z ) = portion of the perturbation potential contributed by sources 
77? ( X, , £ ) = source density on the surface S. 

The function K S (x ; ^, K,t) denotes the potential of an elementary source located at 

(X,^). It is defined as 

s _ -i- V 

M < 1 : K “ ati Joc-x/s))* + W\(r¥, + (z-* ( rs)) 2 / 

(36) 

t s \ M^x-X,(a))-B /(y-y,(s;) a -f O-2,fs))0 
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( 37 ) 


The entire homogeneous solution is given by the sum of these surface integrals as 
(°) v s 

(x, -h U 

Let us now review the behavior of these surface integrals in the limit as the point 
approaches the surface. In particular, it is of interest to examine the resulting velocity dis- 
continuities which appear across the surface. For this purpose we decompose the / integra- 
tion into several parts, separating out for special study a segment of length £ around the 
point on the surface approached by ,2 )as shown in figure 10. It can be shown that all 
discontinuities across the surface arise from the singularities on the small segment £ and that 
the effect of the remaining part of the surface may be disregarded. Also, the effect of surface 
curvature over the small segment E can be disregarded, thereby reducing the problem to that 
of an integral over a planar strip of width E inclined at an angle V to the y -axis. 



After carrying out the integration over the segment with the aid of a Taylor-series 
expansion of m ( X, , J) and ¥(*,,£) about a point on the surface and letting the field 
point (x, y , z ) approach that point on the surface, one finds the following results tor the 
discontinuities across the surface, which are well known. 

, . , , vPP£fiS/D£ Law£J?SIO£ P 

~ = Y’fx^z) - yYx, y /H ;= Jr(x h J)dx, 

— CO 


A u = Y(x, y,i) 

A . 


( 38 ) 
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where 


= discontinuity in the tangential velocity component in the lateral direction 

and 

A ^ = discontinuity in the normal velocity component 

as shown in figure 1 1. 



FIGURE II. -DISCONTINUITIES IN VELOCITY COMPONENTS 

The lateral discontinuities may be decomposed into components in the y and Z 
directions as 


-m(x, y,z)sln V 

a / ( x.y,z) 

s -f 777 ( X, ^ t 2 ) COS (p 


(39) 


It should be noted that in the wake behind a lifting surface only the lateral discontinuities 
arising from the vortex distribution are nonzero. 


11.2 PREFERRED FORM OF THE PARTICULAR SOLUTION 

With these results in mind, we now proceed to write a preferred form of particular solu- 
tion and examine its behavior. The particular solution selected is 
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( 40 ) 


y' (x, y, 2 > = ff ( s} ) r (*• j ^ /x fa 4 ? x, ^ z, i )) d 5 

S 

(K^K 1,*, X,,J))ch5 


+ti s (y-y,(v) ™(x tJ /) & (K 5 (x, y,*, x t/ J))ds 


-If Cx"- x, (s^m(x t J)^(K s (x,^ *, x tJ l)) c/s 


Let us first verify that this is a particular solution of eq. (30). To do this we decompose 
eq. (40) into homogeneous and particular parts. The homogeneous part is 






= y *^s-/£ y, * <c*o/s]+ ^ [/£ X,tK v ds+q x, nK’Js] 


(41) 


The bracketed terms are of the same form as the homogeneous solutions (33) and (35), 
with the vortex and source densities interpreted as (y, ?) , (y, *i) or ( X, Y) t (x, *>)■ They 
are thus homogeneous solutions. The derivatives of these terms appearing in eq. (41) are also 
homogeneous solutions, for it is easily seen by differentiation of eq. (30) that derivatives of 
homogeneous solutions also satisfy the homogeneous equation. Thus the terms in eq. (41) are 
indeed homogeneous solutions. 

The remaining part of eq. (40) can be written as 


H + as -/./*» fps 

- J-fe- [ r&s+tf m k s js\- x §u[ff s t * »<'/5 ] 

- (f*+ f s ) - xfj, (<f v +v s ) 

= y - X 


The latter steps follow from eqs. (33), (35), and (37). It is readily verified by direct substitution 
into eq. (30) that eq. (42) constitutes a particular solution. It thus follows that eq. (40) is also 
a particular solution. 
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Let us now restrict our view to the small strip of width ? , since it contributes all of the 
discontinuous behavior. Define a local coordinate system aligned with the strip, as shown in 
figure 12. 


Z. 



In terms of these coordinates, eq. (40) becomes 

c °* y> // ( tj- - n, ) d ~ dx,^ 

I f V 

-cos yff(x-x, )V d * pffcx-xjr dx, 

t ** f 

* f sir - stn (?- 5 ,) TV 

— Cos ff ( X -x,) m ^rr d x, dr( -h Sin <pfj ( X -X,) nj^r=-dx i 

£ '7 i 


(43) 


where it is understood that y» p now consists only of the contribution from the strip of 
width f . 

It is now a simple matter to deduce the limiting behavior of these terms from a knowledge 
of the basic properties of the integrals which led to the results given by eqs. (38) and (39). For 
example, let us consider the behavior of the first term arising from the discontinuity in . 

It is 

|V£] = Cos // ( V -7/ ) d X ,dJ7/ J + te/»s (44) 
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From eqs. (33) and (38) it follows that 



j 

3 X 





( 45 ) 


This differs from the term appearing in eq. (44) by the factor ( 7 7 - *?i) in the numerator, 
which approaches zero as £ — ► o . Hence the term in eq. (44) is at most ofe) compared 
to ‘fy f X, , and must vanish as C — o , Repeating this type of reasoning for all the 
terms leads to the following conclusions for the velocity discontinuities arising from the 
particular solution: 

F 

Discontinuity in ip x 


Upon differentiating eq. (43) with respect to X and examining the behavior of the 
various terms in the manner described above, one finds only one term that contributes to the 

p 

discontinuity in y; u is 

on- 


(46) 


* - cos f [ _/ 


It follows from eq. (39) and the comment immediately 
nuity across the wake behind the surface is equal to 

[ V* ] " - V [4 J *> ] 


thereafter that the velocity disconti- 
- - [<] 


Thus, the particular solution satisfies the wake boundary condition, eq. (31.) 

P 

Discontinuity in 

p 

The terms from eq. (43) contributing to the discontinuity in *fv are: 
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K J = cos T* d *' d V' j 
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- C^v> ^Iffr^J*.***]-!- l{fa-n,)rT7zr, J X‘ Jt t] 
~[fj^ ^ (x-X/) ^7i, ] J 
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c/x 


-7 


(48) 


r*. ^ ^ 

- Cos f ) (X-K t ) a"^z ^x, 


Discontinuity in —Normal Velocity Component 

" ! r 

The terms from eq. (43) contributing to a discontinuity in <f^ are: 

[%]~ ~ * in V f»T ^ 5 ‘ 3; ) * 2T J 


(49) 
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(49 cont.) 


V/ 

-s<> r[£fr If an*? £r[ ff (?-?,) r%£ 

= - ( X, 7?) Si- >7 </ 

Thus we find a discontinuity in the normal velocity component that is proportional to the 
vorticity distribution of the homogeneous solution. 
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12.0 COMPARISONS WITH KNOWN SOLUTIONS 


To gain confidence in the present theory it is desirable to compare results with some 
known solutions. This can easily be done for the limited class of problems involving planar 
wings in sideslip, since the present theory can then be compared with a first-order expansion 
in 0 of the first-order solution in oi obtained with a skewed planform. Two different prob- 
lems will be considered. 


12.1 INFINITE YAWED WING IN INCOMPRESSIBLE FLOW 

Consider an infinite yawed planar wing lying in the X- ^ plane, with coordinate sys- 
tems and freestream direction as shown in ligure 1 3. The treestream is inclined at the angles 



FIGURE 13.— INFINITE YAWED WING ORIENTATION 


oC and /S to the coordinate planes such that the freestream velocity components in the coor- 
dinate directions are 


V'x ~ ^L. Cos C-os ft 

v = VL COS oc sin 

y ( 50 ) 

V 2 = si rj cc 
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The pressure distribution on this wing can be ascertained by means of simple sweep theory, 
considering the sweep angle to be (Jl + fi ) • This result will then be compared to the pres- 
ent theory, which considers a wing of sweep angle Jl subjected to a sideslip velocity £ 14a • 

To obtain the simple sweep result, we first establish the freestream velocity components 
in the ? , 1} , and S directions. Retaining only first-order terms in aC and £ , a simplifica- 
tion which is consistent with the present theory, these turn out to be 

\T ~ \Z X> CoS -A_ [l ~/3 -to-?? J\- ] 

= \4, co5 -h- T -h- + fi ] (51) 

- **■ 

The pressure distribution is determined by the angle of attack of the normal flow, 
which is 

The latter equality is valid to first order in £ . 

Now the pressure discontinuity across a two-dimensional flat plate of chord c is known 
to be (ref. 3) 



where i/£ and c^o-are the normal freestream component and angle of attack. Substituting 
for and °^zp results in 


P- P~ 

K fio ces 2 J\- [l-/& 4.&>i J 


cos -A. 


+ £ -A- 


] 


HET 

V £ (54) 


which reduces to 

A = y z fv* ^ 1 J ( 55) 

within first-order accuracy. This is the pressure distribution given by simple sweep theory, 
valid to first order in £ . 
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Now let us work out a comparable result by means ot the present theory. The various 
problems, from eq. (28), are 

-h =0 


0(oC) r tp< 


o XX 






3 \ 

wing: - t -a- = - 1 

9?? 


wake 


-Pot 


■■ [Vox ] = ° 

= - 2 Oi <f oK 


(56) 


°W 5 ?,xx * Vyy+Vi 


a? 


= o 


wing: 


=o 

S7 *7 


(57) 


wake: =*• O 


Cp fi = -*/* V?x 

o^/S) ! YW + + = ° 

wing: 4 l — O 

3 n 

wake: ^ ] = - [H* tf] 

* ~ 2 <*(%,* 


(58) 


The total velocity potential for this simple case consists of 

$ * +<*.*+*£•%+* ¥> t + (59) 

The solution to the o(oC) problem is easily derived from simple sweep theory by setting ft 
equal to zero in eq. (55), giving 

a C ?oC = 4 CL Co* -A (60) 

The 0(^5) potential, % , is zero, since all of the boundary conditions are zero. 
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The o(dfi) potential, % , is also zero, since the equation is homogeneous when M -O, 
and all the boundary conditions are zero. Note that, for this particular example, — 0 

in the wake. Hence the only contribution from the o(dfi ) problem comes from the 
expression for C p , which reduces to 


c p«, = - 2 


(61) 


Now, from simple sweep theory, 

tfoy ~ -/i. (62) 


* is related to Cp^ by the expression (eq. 56)) 

C PoC = “ 2 * (63) 

Eliminating <-P 0 x and between eqs. (62), (62), and (63) finally gives 

Cp afi = 'Zaifi <p oX -tan Ji ~ JL Cp a (64) 

Combining eq. (60) and (64) to obtain the total pressure coefficient, one finds that 


ACp = AC Pe< + A 


= 4 o( 


Co : 5 




ACp^ [1 — /3-tanJi 
[ >-£ 


(65) 


Thus, one obtains a term of Ofy) and another of o&fi) • The entire expression is equal 
to the simple sweep theory result of eq. (55), which demonstrates that the present theory 
gives the correct solution for this case. 


12.2 PLANAR WING IN SIDESLIP-SUBSONIC COMPRESSIBLE FLOW 

For this case we shall write the first order in QC solution obtained with a skewed plan- 
form and then expand this solution in powers of the skew angle ft . This will then be com- 
pared with the results of the present theory. 

Consider the wing of planform S shown in figure 14, which lies in the X-y T 
plane. The freestream is directed according to eq. (50). The problem will first be formulated 
in the £ , ^ , 5 coordinate system, in which the freestream has no sideslip component. 
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FIGURE 14.— SKEWED WING PLANFORM 


This problem is 

of*) v 0 -m*) + r S s 

-ZV* = _| 

w,ng: ~Tf 

wake: ~ ° 

where the total potential is 


— a 


« £ [Sf ^5 


( 66 ) 


( 67 ) 


The solution of this boundary value problem can be expressed in terms of an integral 
equation expressing the influence of a distribution of elementary horseshoe vortices. This 
expression is 

a (Sj 7 7 $ ) ~ °< ff*" (It, 7 / J 5 ? * 1 ^5, ( 68 ) 
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where y (% , , % ) is the strength of the bound vorticity and k v is the potential of the 
elementary horseshoe, given by eq. (34) as 





+ 


5 ~ ~5/ 

\lc$- 1,) 1 + (i-m z )(( -n-v.f + 


- (69) 

) 


We now expand this solution in powers of/3 by introducing 


$ ^ X CosJ + y sent ^ x i- fiy 


7? ~ -X3i??/Z +yCos/£ ~ -/3X y ( 70 > 


r- a 


T is expanded as 

ifc.,-n, )**>(*•(*,, 


v,),y. <*>,, v,))*kK <*.&,%)> y,(s,,v.>) m 


Two terms result from the expansion of k y . One js the expression (69) with jf and n 
replaced by X and y , which will be denoted as K ( 


y _ i 


i f- 


X -X, 


4n (y-Vi) 1 *** ( \l(x-k,? + *■ iV 

The other term, of O(fi) , which will be denoted as 3 X 2 , can be written as 

-fx-*.) if] 

Note that this is identical to the form of the particular solution given by eq. (40). 


(72) 


(73) 


With these expansions, the solution (68) takes the torm 


* e) = °< Sf (x, > 2 ,)K*(x, K, V, )<=/*, <=*y> 

-f- Ctfe ^2 ( X, } y i ) K, (X , y > ^ / X, } y,) d X, dyt 

r, ( x, , y,) *2 ( x, y } * > X t , y,) ,<Jy, 
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The boundary conditions become 


wing: f— — — 1 

S 2 


wake 


: [vi] + /*[?*] = ° 


( 75 ) 


Let us now turn to solving the problem by the method presented in this report. The o (oi . ) 
problem is 


(l-M*) 'foxy * £ - 1 - " = ° 




£s 


wing: 


3 77 




(76) 


wake: [> oX ] = ° 

The 0^4) solution is identically zero, because of the zero boundary conditions. The o(atfi) 
problem, however, is 


*• x 

wing: ^ ^ — q 

3 h 


(77) 


wake: = - |>„ y ] 

The total velocity potential is 

k~ ^ r * + fi i + 2 + * ^ ^ ] (78) 

The solutions to these problems are 

(* tJ yt) K ?( X >?J i , y ' > yi) J y , ( ?9 > 


$$ *4 (x • , y>) X > y > X f > it) d y, (80 ) 

+ ff s r 2 (*,, i,) K% (X > **> y,)dx, Jyj 

where *3 and *4 are established by the boundary conditions. The latter term is composed 
of a particular solution involving J< * (see eq. (73)), and a homogeneous solution containing 

<r • 
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Having obtained a solution by both methods, it remains to compare them. First note 
that the solution ( 74 ) is composed of an o(o( ) term and two o(<* fi) terms, which would be 
identical to the solution by the present method, the sum of eqs (79) and (80), provided that 
it can be shown that 


t , - ^3 

(81) 


(82) 


This can be done as follows. Let us write the term o< as 

<*</«• « ft* + <*/* ft* 

C( 

where oi y/* is equal to the first right-hand term of eq. (74) and cC& ^ equal to the 
remaining terms. If we can show that S\-^3 and , then / 0 and 

will be equal to their counterparts ft and ftt obtained by the present method. Substi- 
tuting the expression (82) into the boundary conditions (75) gives 


wing: 

3 ftf 

+ 

4 

» ft* 

— - 1 

3 1 

3 z 

wake : 


i + 



1 +/[<] + 


-A o(F)~0 


(83) 


Equating orders of fi gives 

c( 

O(oC) • wing: ^ ~ * 

e> t 

: &*] * 


wake: 


O (elfi )*• win S- ~ ° 


wake: 




(84) 


(85) 


In comparing these with the boundary conditions given by eqs. (76) and (77), it is 
obvious that the boundary conditions determining 2 T { and £ 2 . are identical to those 
determining and 9 a , and hence that 1 — £3 and . Thus the equiva- 

lence of the solutions by the two methods has been proved. 
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13.0 CONCLUDING REMARKS 


A theory has been developed for predicting the aerodynamic properties of an airplane 
in sideslip. The airplane geometry was defined in terms of small parameters governing the 
angle of attack, angle of sideslip, wing camber, wing thickness, body radius variation, and 
body camber. These small parameters formed the basis for a perturbation expansion approach 
wherein the velocity potential and boundary conditions were expanded as asymptotic series 
in powers of the small parameters. All tirst-order terms and those second-order terms involv- 
ing the interaction between sideslip and the other small parameters were retained. 

The development produced a set of first-order homogeneous boundary value problems 
giving the influence of each small parameter, and a set of second-order inhomogeneous prob- 
lems giving the interactions between sideslip and angle of attack, sideslip and camber, etc. 

It was found that certain of the inhomogeneous terms appearing in the second-order problems 
were always of higher order than those producing the dominant effects and hence could be 
deleted. This led to a simplification of the second-order problems such that generalized par- 
ticular solutions could be obtained. A particular integral was found which has the property 
of automatically satisfying the second-order boundary conditions on all trailing vortex sheets. 
With this, the remaining second-order homogeneous part of the solution can be constructed 
by the same methods used to solve the first-order problems. The theory was checked by com- 
paring with known solutions for an infinite yawed wing and for a skewed flat wing. 

A major conclusion is that the dominant effects influencing the aerodynamics of side- 
slip can appear either as first-order or as second-order terms, depending primarily on the 
magnitude of the dihedral of the lifting surfaces. For configurations having large amounts of 
dihedral, the first-order terms dominate. Sideslip aerodynamics will then depend primarily on 
the angle of sideslip alone, and will not be strongly influenced by interactions between side- 
slip and angle of attack or between sideslip and the other parameters. For such configurations, 
the aerodynamic description given by conventional first-order theory is sufficient. 

However, if the dihedral is small, then first-order theory fails to give the dominant effects. 
For such configurations, the second-order terms involving the interaction between sideslip and 
angle of attack, etc., are major contributors to sideslip aerodynamics, and it is essential that 
they be taken into account. 

Another outcome is that highly accurate numerical solution procedures are necessary in 
any calculation of second-order sideslip aerodynamics. This requirement arises from the singu- 
lar nature of the kernel of the particular integrals for the second-order problems. The kernel 
is more singular than its homogeneous counterpart, and integration by parts can be used to 
show that the numerical value of the particular integral depends strongly on the local gradients 
of the coefficients involving the first-order source and vortex strengths. This behavior can also 
be seen in the alternate form for the particular integral, 

v ,p = y f x re?; - * 
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from whence it is obvious that an accurate knowledge of the velocity derivatives of the first- 
order solution is required for the evaluation of aerodynamic forces and pressures contributed 
by the second-order particular solution. It is apparent, then, that the present theory must be 
accompanied by a numerical solution procedure that meets requirements more stringent than 
those ordinarily needed in a first-order analysis. For an ordinary first-order analysis the aero- 
dynamicist is usually satisfied with a numerical solution procedure that provides accurate 
forces and pressures (velocities). In contrast, the present second-order theory must be accom- 
panied by a numerical method capable of producing accurate velocity gradients in the first- 
order solution. 


The Boeing Company 
P. O. Box 3707 

Seattle, Washington 98124 
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